How to make a bilayer exciton condensate flow 
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Bose condensation is responsible for many of the most spectacular effects in physics because it 
can promote quantum behavior from the microscopic to the macroscopic world. Bose condensates 
can be distinguished by the condensing object; electron-electron Cooper-pairs are responsible for 
superconductivity, Helium atoms for superfluidity, and ultracold alkali atoms in vapors for coherent 
matter waves. Electron-hole pair (exciton) condensation has maintained special interest because it 
has been difficult to realize experimentally, and because exciton phase coherence is never[l| perfectly 
spontaneous. Although ideal condensates can support[2| an exciton supercurrent, it has not been 
clear[3| how such a current could be induced or detected, or how its experimental manifestation 
would be altered by the phase-fixing exciton creation and annihilation processes which are inevitably 
present. In this article we explain how to induce an exciton supercurrent in separately contacted 
bilayer condensates, and predict electrical effects which enable unambiguous detection. 

PACS numbers: 



The order parameter of an exciton condensate is 
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\^(r)\exp(i(j){f)) = (ipl(r)iph(r) ) = p(h,f;e,r) 

(1) 

where i/j> and i/j are electron creation and annihilation 
operators, (f>(r) is the condensate phase, the labels e (elec- 
tron) and h (hole) refer to the states between which phase 
coherence is established (nearly!) spontaneously, and 
p(h, r, e, r) is the anomalous density matrix. Microscopic 
considerations suggest [4| that spontaneous coherence is 
likely between a conduction band with occupied states 
inside a Fermi surface and a valence band with occupied 
states outside a nearly[5| identical Fermi surface. Part of 
the reason that exciton condensation has not been easy 
to realize is that sufficiently perfect nesting between con- 
duction and valence bands is unlikely to occur naturally. 
The systems of interest here are artificially fabricated 
bilayers in which the electrons and holes are in well sep- 
arated two-dimensional electron systems(2DESs), either 
semiconductor quantum wells 0, 0] or graphene layers @ 
separated by a dielectric barrier. The dielectric barrier 
reduces the strength of exciton creation and annihilation 
processes, and gate control of the density in each layer al- 
lows electron and hole band Fermi surfaces to be tuned to 
the same area. Although simple to describe, this quan- 
tum engineering is difficult [1, @, to execute success- 
fully. Exciton condensation has so far been realized [ici| 
only in the quantum Hall regime in which band disper- 
sion is irrelevant allowing spontaneous coherence to occur 
between spatially separated conduction bands, or spa- 
tially separated valence bands, under circumstances that 
are achieved routinely. The considerations explained in 
this article apply to quantum Hall exciton condensates 
in the Corbino geometry [Til [T3|. in which current flows 
across the 2DES bulk, but not directly to the Hall bar 
geometry[l3j in which currents flows along the 2DES 
edge. 

In their pioneering work on exciton Bose condensation 
Blatt et aZ.[!| argued that because an exciton is neu- 
tral, condensation can not lead to spectacular electrical 



effects. Experimental studies of quantum Hall exciton 
condensates have already made it clear [l(| E3] that this 
pessimistic conclusion is not valid. The key technical 
capability not anticipated in 1962 is the possibility of 
making independent electrical contact [HJ to the electron 
and hole parts of the condensate. In this article we ex- 
plain how condensation leads to a reorganization of the 
low-energy charged fermion degrees of freedom which is 
responsible for dramatic changes in the transport prop- 
erties of separately contacted condensed bilayers. 

Spectacular electrical effects in exciton condensates are 
enabled by the possibility of exciton supcrflow. The key 
issue which arises in addressing these phenomena theo- 
retically is understanding how a supercurrent of neutral 
excitons can be driven by electrochemical potential dif- 
ferences. We consider a bilayer with contacts to both left 
(L) and right (R) ends of the separate quantum wells (top 
(T=e) and bottom (B=h)) between which coherence is 
established, as illustrated schematically in Fig.[[l]. (In 
the case of the Corbino geometry transport of quantum 
Hall condensates L and R refer to the outer and inner 
edges of an annular 2DES.) For the sake of definiteness 
we focus our attention on voltage biased transport; our 
analysis is easily extended to the case of current biased 
transport and to systems with additional voltage probes 
along the sample length. The observables in this trans- 
port geometry are the currents and voltages in the leads. 
Exciton condensation alters transport properties by in- 
troducing an anomalous symmetry-breaking field in the 
fermion quasiparticle Hamiltonian (more technically in 
the one-particle Greens function) which enforces inter- 
layer phase coherence, opens up a gap at the Fermi level, 
and allows charge to move freely between layers. In 
mean-ficld-thcory, the quasiparticle amplitude for tun- 
neling between top (T) and bottom (B) layers is: 

(T,r'\H HF \B,r) = V C (\P -f\) p{T,?;B,r) (2) 



where V c (r) = e 2 /e\/r 2 + d 2 is the inter-layer Coulomb 
interaction and p(T,r';B,r) is the inter-layer compo- 
nent of the density matrix. (p(T,r';B,r) is non-zero 







FIG. 1: Schematic illustration of a separately con- 
tacted bilayer exciton condensate. L and R represenl 
the left and right ends of a quantum well system or the outei 
and inner edges of an annular quantum Hall bilayer conden- 
sate. Currents through the four leads can be tuned externallj 
by varying the three independent voltage differences. Oui 
convention for positive currents is indicated by the arrows at- 
tached to the leads. Transport properties are simplest wher 
there is a gap (indicated by white fill in this schematic illus- 
tration) in the bulk. 
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FIG. 2: Loop and load geometries. Fig. a and b represent 
respectively loop and load geometries. For bilayer exciton 
condensates only the load geometry b permits steady state 
exciton superflow. 



only in the condensed state.) These quasiparticle tun- 
neling terms are present even when the microscopic elec- 
tron Hamiltonian does not include inter-layer transfei 
processes. In a mean-field theory microscopic particle- 
number conservation in each layer is recovered when 
the exciton pair potential is calculated self-consistently; 
inter-layer quasiparticle currents arc then precisely com- 
pensated by exciton condensate supercurrents. Steady 
state exciton supercurrents are possible, but only when 
the lead voltages are chosen to satisfy a condition that 
we explain below. Our main predictions are qualita- 
tive in nature. The quantitative properties of a partic- 
ular exciton condensate system will depend on the de- 
tails of its Hamiltonian, including its disorder potentials. 
We illustrate some of our ideas on transport in exci- 
ton condensates below by considering a one-dimensional 
toy model with short-range inter-layer interactions, us- 
ing the NEGF[n| method to evaluate steady-state non- 
equilibrium density-matrices. 

When p(T,f; B,r) has a spatial phase gradient, the 
quasiparticle Hamiltonian supports [2| non-zero currents 
even when the system is in equilibrium, i. e. when all four 
chemical potentials are identical. Because the exciton 
condensate cannot carry charge, the condensate current 
(CC) has onlyQ a counterflow component, i.e. 

3? C (x)+3% C (x)=0. (3) 

The total current in each layer can be expressed as the 
sum of the condensate current and a quasiparticle cur- 
rent (QC) driven by lead voltage differences. The QC 
conduction is related via the Landauer-Buttiker formula 
to the transmission coefficients of quasiparticle waves in- 
cident from the various leads. In order to drive steady 
state condensate current, the gate voltages must be cho- 
sen so that the counter-flow component of the QC is not 



spatially constant: 
ld x (j? C (x) - f B C {x))] = -[d x {§ C {x) - jg C (*))] (4) 

Condensate currents are induced by space-dependent 
quasiparticle counterflow currents. This relationship be- 
tween quasiparticle (QC) and condensate (CC) counter- 
flow currents follows from the separate conservation of 
charge in each layer. Since condensate currents cannot 
enter the leads, they must be present in the bulk when 
across the sampleirL + Ibl ^ or Itr + Ibr ^ 0. (Here 
we use the conventions of Fig.[[T] for the lead current 
signs.) For every electron that is reflected from the top 
layer to the bottom layer at the left edge of the sample, 
an exciton must be launched. It follows that (Itl + Ibl) 
is the exciton supercurrent emitted from the left side of 
the sample and that (Itr + Ibr) is the exciton super- 
current absorbed on the right. Steady state condensate 
currents are possible only when these two quantities are 
equal: 

Itl + Ibl = Itr + Ibr- (5) 

In the linear response regime, the four lead currents are 
proportional to the three independent electrochemical 
potential differences. Eq.(|4]) places a restriction on the 
two independent difference ratios. 

One convenient way to reduce the lead-voltage-space 
dimension experimentally is to connect two contacts with 
a resistor R, for example in either the loop geometry or 
the load geometry illustrated in Fig. ([5]). In the load ge- 
ometry, the contact voltages arc limited to the surface on 
which Ibl = Ibr = (Vbl — Vbr)/R- The consequences 
of this resistive link between contacts are easily antici- 
pated in the most common limit in which the gap due to 
condensation is large enough or disorder strong enough 
to prevent quasiparticle current conduction across the 
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sample. Bulk quasiparticle currents are expected to be 
negligible for Corbino geometry transport in the quan- 
tum Hall regime, and we anticipate that they will also be 
negligible in zero-field bilayer exciton condensates when 
these are ultimately realized. When quasiparticles can- 
not flow across the sample, Itl — Ibl and Itr = Ibr- 
It follows that Eq. ( O is always satisfied in the load ge- 
ometry. Because it is possible to induce an exciton su- 
percurrent in the bulk of the bilayer, a large charge cur- 
rent can flow through the circuit with a resistance due 
only to the load resistor and contact resistances at the 
two ends of the sample. In the case of quantum Hall 
bilayers, we predict a Corbino resistance for the load ge- 
ometry which remains constant as temperature T — > 
on bilayer quantum Hall plateaus which are due to exci- 
ton condensation [14| . This behavior is in stark contrast 
with the dramatically increasing resistance expected at 
low temperatures on plateaus not associated with exciton 
condensation. 

Although the loop geometry of Fig. ( [J) is naively com- 
patible with steady state counterflow currents, we predict 
that it cannot support steady state exciton superflow be- 
cause it does not guarantee that the the exciton super- 
current emitted at the left end of the sample is identical 
to that absorbed at the right end. Indeed for Itr> the 
loop geometry will lead to Vtr > Vbr and therefore to 
quasiparticle current flowing in the wrong direction from 
top to bottom. 

In the load geometry, the effective two-probe quasipar- 
ticle conductances at left and right Gl,r — (e 2 /h)TL,R 
depend on the details of these contacts and on the quasi- 
particle Hamiltonian near the sample ends. Taking the 
TR lead as ground, the voltages on the other leads in 
this circumstance are V T l = I (PV 1 + T R ^(h/e 2 ) + R), 
Vbl = I {R+T R l (h/e 2 )), and V BR = I (h/e 2 )^ 1 where 
Itl = Ibl = Ibr = Itr = I- When the load resistor 
R is small, the current which flows through the bilayer 
systems is limited only by the quantum contact resis- 
tances (h/e 2 )T£ R . Even when the microscopic electrons 
tunneling amplitude is zero, quasiparticles move freely 
between layers allowing charge to be conducted across 
the highly resistive bulk by exploiting the parallel load 
resistor channel. 

In Fig.( [3]) we plot numerical results for the self- 
consistently calculated quasiparticle and condensate cur- 
rent distributions in a one-dimensional tight-binding- 
model bilayer system with an electron band in the top 
layer and an equal density hole band in the bottom layer. 
The inter-layer tunneling amplitude in the quasiparticle 
Hamiltonian of this model is 



(T, I'\H HF \B, I) = Srj V p(T, I'; B, I) 



(0) 



where V is the short-range inter-layer interaction 
strength and /, I' is the site indices. Because of the 
energy gap induced in the quasiparticle Hamiltonian by 
this tunneling term, the quasiparticle current is deflected 
to flow between T and B leads at both left and right. 
The numerical results in Fig.( [3]) have been obtained 
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FIG. 3: ID Toy model self-consistent current distribu- 
tions. Sites 1 through 300 label the top layer left to right 
while sites 301 through 600 labels the bottom layer. The 
site separation is 0.08 Bohr radii. The black solid lines, red 
dashed lines, and blue dotted lines represent the total current, 
the condensate current (CC) and the quasiparticle current re- 
spectively. In a,VrL = 1.0, Vtr = —0.5, Vbl = Vbr = 0.0 
and the bare interlayer tunneling amplitude A t = 0.0. (Volt- 
ages are in Rydberg/e) The total (CC+QC) current is con- 
stant and the QP current is very small near the system center. 
h,V TL = 0.5, Vtr = -0.5, Vbl = 0.0, Vbr = 0.0, A t = 0.005. 
When bare tunneling occurs the total current is no longer 
a constant, but the transport characteristics in loop geome- 
try do not change in a qualitative way. c,Vtl = 0.5, Vtr = 
-0.5, Vbl = 0.0, Vbr = -0.1, A* = 0.005. Steady state so- 
lutions with net exciton creation rates are possible because 
of bare electron interlayer tunneling. The unit of voltage is 
Rydberg/ e. 



by using mean-field theory (Eq.( \§§) to calculate the 
quasiparticle Hamiltonian given the density matrix, and 
the NEGF formalism[l6[ to evaluate the non-equilibrium 
density matrices given the lead voltages and the quasi- 
particle Hamiltonian. In the top panel of Fig.( [3]) 
,Vb(l,r) = ((Vtl + V T r) ± V')/2. Because the left 
and right sides of our model system are identical this 
choice produces Ibl = Ibr, corresponding to the load 
geometry. (V is related to the load wire resistance by 
V'/{V T l - V TR ) = R/(R + TV/ 1 + T^ 1 ).) When the 
non-equilibrium mean-field calculation is carried to self- 
consistency, the phase of the order parameter develops 
a spatial gradient and the system carries a steady state 
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condensate current. 

The total current is constant in each layer as required 
by charge conservation. In the ordered state, the low- 
energy charged degrees of freedom are reorganized from 
bare electrons which cannot transfer between layers to 
exciton-condensate quasiparticles which occupy the two- 
layers simultaneously and coherently. This profound re- 
organization of low-lying charged excitations decreases 
the resistance exponentially compared to the case in 
which the gap in the quasiparticle spectrum is not due 
to exciton condensation. In the quantum Hall case for 
example, we predict load geometry resistances that are 
orders of magnitude smaller at v — 1 (the exciton con- 
densation case) than at v = 2. For the loop geome- 
try, no self-consistent steady-state solutions of the non- 
equilibrium self-consistent field equations exist. In this 
case we anticipate enhanced noise due to order-parameter 
time dependence and dramatically higher steady state 
time-averaged resistance. 

As emphasized in a fundamental early paper by Kohn 
and Sherrington p|, phase coherence between different 
bands in a solid can never be completely spontaneous. 
In the case of bilayer condensates, inter-layer electron 
tunneling, which creates or annihilates excitons, is ex- 
pected to be the dominant process which fixes a preferred 
phase and violates separate charge conservation in each 
layer. As illustrated in panel b in Fig.( [3]), however, 
weak inter-layer electron tunneling has little effect on 
transport properties. Most experimental anomalies as- 
sociated with bilayer exciton condensation require only 
that the quasiparticle tunneling amplitude be dramati- 
cally enhanced compared to its bare values. As we now 
explain the main consequence of bare-electron inter-layer 
tunneling is partial relaxation of the current-conservation 
condition, Eq.([5]). 

The role of a bare electron inter-layer tunneling am- 
plitude is most simply discussed using the minimal field 
theory model[2| of a bilayer exciton condensate: 



E 



dr 



2m* 



|V0| S 



M At cos 9 



(7) 



where 9 is the inter-layer coherence phase, A t is the bare- 
electron tunneling amplitude, ,to* is the effective mass 
in the Landau- Ginzburg type theory, M = 2\p(T,B)\ is 
the condensate order parameter, and p s is its superfluid 



density. In this description j = (Hp s d x 9)/m* is the ex- 
citon condensate supercurrent. Minimizing this energy 
functional leads to a sine-Gordon equation from which it 
follows that 



\\ 2 (d x 9) 2 +cos0 



0. 



(8) 



where A 2 = (p s /M)(h 2 /to* A t ) is the model's Josephson 
length. From the constant-of-motion in square brack- 
ets in Eq. ( [8]) , we obtain an explicitly expression for the 
difference between the condensates currents at opposite 
ends of the sample: 



3% 



jI 



x 2 fhp s 
(cos9 L - cos 6^) — — 

\ z \ TO* 



(9) 



Defining J = (]r + j'l)/2 and Sj = jr — jl, we conclude 
that steady state collective currents are possible provided 
that the exciton currents injected at left and absorbed at 
right differ by less than 



< A, 



2p s M 1 

TO* J 



(10) 



The limit on Sj is closely analogous to the critical cur- 
rent predicted[l3| for inter-layer tunneling. For contact 
resistances are ~ h/e 2 we predict that large conductances 
will occur when Vbr is varied by less than ~ (h/e) Sj max 
from its load resistance value. In our one-dimensional toy 
model for example (j L = Itl/z, 3r = Itr/z), Eq.([10} 
implies that Sj max ~ 0.1(Rydberg//i) 

for A t = 0.005 (Rydberg) and total current j = 
0.5(Rydberg//i). Indeed, our NEGF calculations find 
a self-consistent solution, illustrated in Fig.( [3] c), for 
the voltage configuration Vtl = 0.5, Vbl = 0.0, Vtr = 
— 0.5, Vbr = —0.1 (in (Rydberg/e) units) (corresponding 
to Sj ~ 0.1(Rydberg/ft,)) but not for larger Sj in quali- 
tative agreement with this prediction. 
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